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Abstract 

We consider a sequence of linear hyper-elastic, inhomogeneous and 
fully anisotropic bodies in a reference configuration occupying a cylindri¬ 
cal region of height e. We then study, by means of F-convergence, the 
asymptotic behavior as e goes to zero of the sequence of complementary 
energies. The limit functional is then identified as a dual problem for a 
two-dimensional plate. Our approach gives a direct characterization of 
the convergence of the equilibrating stress fields. 
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1 Introduction 

The equilibrium problem for a linear hyper-elastic body may be suitably studied by 
means of several variational formulations, like the principle of the minimum potential 
energy (primal formulation) and the principle of minimum complementary energy (dual 
formulation). In the former formulation the unknown is the displacement vector field, 
while in the latter the stress tensor field is to be found. Other variational formulations, 
the so called mixed formulations, take simultaneously as unknowns the displacement 
and the stress vector fields, see for instance [I]. 

In the last three decades, starting with the work of Ciarlet and Destuynder [2], 
these variational problems, or their extremality equations, have been used, in con¬ 
junction to some asymptotic techniques, to justify/derive models for thin structures 
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starting from the three-dimensional theory. At the early stages of this prolific line of 
research mixed formulations were adopted, while after the asymptotic techniques have 
been refined the research have been focused almost exclusively on the study of some 
form of the primal formulation. Within this line of research, the Kirchhoff-Love theory 
for homogeneous and isotropic plates has been justified by means of T-convergence by 
Anzellotti et al. [5] and by Bourqin et al. 0. These results have been generalized in 
several directions: for linear plates with residual stress 00 , for elasto-plastic plates 
000 , Reissner-Mindlin plates and non-linearly elastic plates EHU. 

Respect to the existing literature a different route has been taken by Bessoud 
et al. in | 15 |. These authors consider a system of two elastic materials glued by a 
thin and strong material between them and by means of the complementary energy 
they study the asymptotic behavior of the system of materials as the thickness of the 
gluing material goes to zero. In the limit problem a material surface, endowed with 
an appropriate elastic energy, replaces the thin layer. 

We here consider a sequence of linear hyper-elastic, inhomogeneous, and fully 
anisotropic bodies in a reference configuration occupying a cylindrical region of height 
e. We then study, by means of T-convergence, the asymptotic behavior as e goes to 
zero of the sequence of complementary energies. The limit functional is then identified 
as a dual problem for a two-dimensional plate. 

While variational limits of primal problems characterize the asymptotic behavior of 
the minimizing displacements, the study of the asymptotic behavior of the complemen¬ 
tary energies characterizes the convergence of the equilibrating stress fields. Besides 
the use of this novel approach for the deduction of plate theory, our work deals with 
fully anisotropic and inhomogeneous materials, case that has not been studied in this 
full generality before, not even by means of the primal formulation. This kind of 
generality on the constitutive equations has been used to derive linearly elastic beam 
theories, within the primal formulation framework, for instance in ph] il71 [T 8 l 119) . 

The paper is organized as follows. In Section 0 we review some function spaces 
that will be useful in the rest of the paper, while in Section 0 the primal and dual 
formulation of the problem considered are stated. The dimension reduction problems 
are classically rescaled on a fixed domain, this is done in Section 0 In Section 0 the 
Gamma-convergence analysis is carried on, and in Section0the obtained Gamma-limit 
is written on a two-dimensional domain. 


2 Preliminaries 

Let 0 C R 3 be an open, bounded set with Lipschitz boundary dfi, and let T be an 
open subset of aeQ We denote by 

H 1/2 (T- R 3 ) := {v : 3u € tf^^R 3 ) s.t. = v on T}, 

where 7 : H 1 (Q\ R 3 ) —¥ R 3 ) denotes the trace operator, and we equip it with 

the norm 


IMIffi/2(r) : = infllMItfHfi) : u € H 1 ^; R 3 ) and 7 u = v on T}. 
The dual of tf 1 / 2 (T; R 3 ) shall be denoted by tf- 1 / 2 (T; R 3 ). We let 

^ 00 2 (T; R 3 ) := e H 1 / 2 (T;R 3 ) : v £ H 1/2 (dQ; R 3 )}, 

1 Note that no regularity assumption is made on the open set T. 
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where v denotes the extension by 0 of v to and we equip it with the norm 

IMI w i/ 2 ( r) := IMIiH/^an)- 

The spaces Id 1 / 2 (r;R 3 ) and H^ 2 {T- R 3 ) are delicate spaces, see e.g., see [20] Chapter 
I, §11 and 12, and [21] Chapter 1, §1.3.2. Note that the space denoted here, and in 
[ 20 ] , by Hqq 2 (T) is denoted by VF 2 1 ' / 2 (r) in [ 21 ] . 

Thanks to these spaces, for a distribution / £ H^ 1 ^ 2 [dVL\ R 3 ) defined in the whole 
boundary we may define its restriction to T, denoted by /|r £ (H^ 2 (T;S. 3 ))', in 
the following way 

(f\ r ’ v ) h^ 2 (t) := (f’V) H i/2 on) for ever y v e #dd 2 ( r ; R3 )- 

The space 

H(div,n) := {T £ L 2 (fl;R 3 y x ^) : divT £ T 2 (n ; R 3 )}, 
equipped with the norm 

l|2"lllf(div,r 2 ) : = l|r|| 2 2 ( n ) + ||divT|| 2 2 (Q), 

is a Hilbert space. It is well known that there exists a continuous linear mapping 
7 n : H(div, H) —> H~ 1 ^ 2 (d 0;R 3 ) such that 

/ T ■ Vudx = — divT • uda: + ( 7 „T, 7 u) if i/ 2 (a m, (1) 

in in 

for every T £ H( div, H) and u £ (f2; R 3 ). Hereafter we shall simply write Tn in 

place of 7 n T. 

From ©it follows that for T £ H( div, fl) and for every u £ H 1 (f2; R 3 ) with yu = 0 
in H 1/2 (dQ \ F; R 3 ) we have 


/ T ■ S7udx = — divT • udx + (Tn|r, 7 «) 
in in 


r od 2 ( r )’ 


( 2 ) 


since 7 M £ HgQ 2 (r;R 3 ). 

Hereafter, if no confusion shall arise, we shall simply write Tn also for the restric¬ 
tion Tn|r and we shall drop the use of 7 to denote the trace, i.e., we shall write u for 
the trace -yu. 


3 The unsealed problems 

Let u> be an open bounded domain of R 2 with Lipschitz boundary duj, and for e £ (0,1], 
we set 

Q e := w x (-e/ 2 , e/ 2 ). 

Let 8du) and c/vtu be unions of finite numbers of open connected subsets of dui such 
that 

Odoo (~l di vlo = 0, 8doo U 8nu = dcu, 8 du ^ 0. 

We set 

do^ e ■= 9du x (—e/ 2 , e/ 2 ), diyQ e := \ dijfT. 


3 





We consider Q £ as the region occupied by a linear hyper-elastic body in the reference 
configuration. Let C e £ L°°(Q, £ \ R 3 j?m X3x3 ) be the elasticity tensor, which we assume 
to be coercive, of the elastic body considered. By writing C e £ Z/°°(fi E ; R 3 ^ x3x3 ) we 
mean that 

Cf jki = = Cfji/b- 

The sets 8dTL £ and <9jvf2 E are the parts of the boundary of fl E where Dirichlet and Neu¬ 
mann boundary conditions are imposed, and we denot(0 by f £ £ (Hqq 2 (dN^ e R 3 ))' 
the surface loads, and by g £ £ H 1 ^ 2 (dD^i e ', R 3 ) the imposed displacement on 9 d 12 e . 
Since every function in H 1,/2 (<9£>fi e ;R 3 ) is the trace of a function in Lf 1 (f2 e ;R 3 ), we 
also denote by g £ this latter function. 

We further denote by b £ £ L 2 (f2 E ;R 3 ) the body forces. 

Remark 3.1 Let u)± be the upper and lower bases of the cylinder fl E and let be 

the Neumann part of the lateral boundary, i.e., 

lo± :=w x {±e/2}, and OnTiI := <9jv<u x (—e/2, +e/2), 


so that U u)L U OnQ\ is di vf2 E up to a set of zero two-dimensional measure. Let 
f± £ H~ 1 ^ 2 (u}±', R 3 ), and ff £ R 3 ). Then, f £ define (0, for every v £ 

Htf(d jvn e ;R 3 ), by 


+ (/^ i '3|9jvf*|)ir 1 / 2 (s N -n|)i 


(3) 


is an example of a force that can be used as f £ . Moreover, in the case that f± £ 
Z/ 2 (a>±;R 3 ) and ff £ L 2 (9jvO|; R 3 ), the duality in ((3| is nothing but the sum of three 
integrals. Note however that there are forces in (Hqq 2 (dN^l e R 3 )) / that are more gen¬ 
eral than f £ defined by ©■ 


Remark 3.2 Given f £ £ (RgQ 2 (0ivfi £ ; R 3 )) 7 and b £ £ L 2 (f2 E ;R 3 ), as above, there ex¬ 
ists (looking e.g. for G £ as the symmetric part of the gradient of an unknown function) 
a G £ £ H( div, fi e ) such that 

( divG £ + b £ = 0 in L 2 (Q e - R 3 ), 

\ G £ n = f in(H^ 2 (d N Q £ )Y. 

Then, from ©. the work done by the loads can be simply rewritten as 

jj- vdi + </*,«>„„.,= £ c* . Evdi, (5) 

for all v G ff 1 (f2 e ;R 3 ) such that v = 0 on <9d^ £ . 

Also the converse is true: given G £ G div, Q e ) there exist f £ G (Hqq 2 M 3 )) 7 
and b £ G L 2 (Q £ ;M 3 ), defined by (|4|). for which (|5j holds. 

2 Throughout the paper the notation • refers to quantities defined on or on parts of its 
boundary. 

3 Throughout the paper the notation 7 refers to quantities which are examples of the general 

case. 
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Therefore the description of the applied loads may be done indifferently either by 
means of the body and surface forces b e and f e , or by means of the tensor field G e . 
Both approaches present some advantages and some disadvantages. For instance, it is 
not necessary to assume G e £ H(div, fl e ) but it is enough to have G e £ L 2 (fl e ; R 3 .^). 
We consider hereafter both representations simultaneously. In Remark below, we 
explain why we consider both type of forces. 

We consider, in the spirit of Remark 13.21 also “generalized forces” described by a 
tensor field H e £ L 2 (fT; R 3 *^). 

The problem of linear elasticity can be written as: 

' w e £ H\n e -R 3 ),w £ = g e in H 1/2 (d D n e -, R 3 ), 

, J^C e Ew e -E^dx = j^H e -E^ + b e -Tpdx + if'J)^^^, (6) 

for every %/j £ A e , 

where Eut e denotes the symmetric part of the gradient of w e , and A e the set of 
admissible displacements defined by 

A e := {6 £ H\ f2 s ;R 3 ) : v = 0 on d D n e }. 

Since g s simultaneously denotes a function in H 1 (Q, £ ;R 3 ) and its trace, ([6]) can be 
rewritten as: 

' u e :=«f -g e £ A e , 

< J^C e Eu e -E$dx = f 'F* (7) 

for every ip £ A e , 

where, for notational simplicity, we denote 

F e :=H e + C e Eg s . (8) 

As it is well known, the solution u e of 0 may also be found by minimizing the 
total energy F E : A e —¥ R defined by 

C e Ev-Evdx-J^ F e -Ev + b e -vdx-if^v)^^^, 

that is 

P{u e ) = inf P(v). 

v£A £ 

This variational problem is called Primal Problem. 

We now introduce the dual problem. 

From the inequality 

0 < t e {E - (C £ ) _1 ,$) ■ (E - = C e E -E-2E-S + (C e ) _1 S • S, 

which holds for every E, S £ R 3 ^, it follows that for every v £ H 1 ( R 3 ) 

-f C £ Ev ■ Ev dx = max [ S ■ Ev - S ■ Sdx 
2 Jne sen 2 (n e ;M^) Jne 2 
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with the max achieved for S = (C £ Ev. Thus we can rewrite the direct problem as 

inf max C £ (v,S), 
veA? sen 2 (n £ ;Ksym) 

where the Lagrangian £ £ (v, S) is defined by 

C £ (v, S) i(C s )”^ ■ S — F e ■ Ev -b £ -vdx — (f £ ,v) H i/^^. 

Since £ e satisfies the assumptions of the min-max Theorem (see e.g. |22 | p. 176 Propo¬ 
sition 2.4 and Remark 2.4), it follows that 

inf max C £ (v,S) = max inf C e {v,S) 
veA? seL 2 (seL 2 (a^-,mif^)veA s 

= max ( f — i(C E ) _1 5 • Sdx 

S€-L 2 (Q e ;Ksym) 4 ns 2 

+ inf f {S-F £ )-Ev-b £ -vdx-{f £ ,v) 1/2 ), 

veA? Jn e 00 { N ’ 

but, from m we deduce that 

inf [ (S-F £ )-Ev-b e -vdx-(f e ,v) Hl/2 
veA £ Jqb h oo Wn** ) 

= inf ( [ (-div{S-F e )-b e )-vdx + {(S-F e )n-f,v) 1/2 ) 

J 0 if S £ S £ , 

\ — oo otherwise, 

where we denote by 

S £ := {S' € i 2 (n s ;R^) : di v(S-F £ ) + b £ = 0 in L 2 {Q. £ - R 3 ) and 

(S - F e )n - f £ = 0 in (H^ 2 (d N Q £ ; R 3 ))'} 

the set of admissible stresses. 

By defining the dual energy by 

F e *(S) := i [ (C £ )~ 1 S-Sdx, 

2 in® 

it follows that 

P(u e ) = inf P{v) = - min F £ *{S) =: -F £ *(T £ ), (9) 

fie^fs sess 

and that T e = C £ Eu £ . 

The minimization problem 

min F £ * (S). 

se5 e 

is called Dual Problem. 

Remark 3.3 Note that the stress d £ := C £ Ew £ associated to the solution ut £ of ©, 
see also ©, is given by 

a £ = T £ +C £ Eg £ . 
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4 Rescaled problems 

We now rescale the problems introduced in Section [3] to a domain independent of e. 
To this end, we set 


:= f2i, 8d LI := 8nLI i , 8dLI ’■= 8nLI i . 
We define the change of variables p £ : SI —> fl e by 

P £ (X 1,X2,X 3 ) ■■= (X!,X2,eX3), 


and we let 

P e : = Vp e = diag (1, l,e). 
For v : Q, e —> R 3 we define v : —> R 3 by 


P E £ 

V o p , 


( 10 ) 


so that 

Vv = P £ (Vv ) op £ P £ , and Ev = P £ (Ev) o p £ P £ . 

We denote by 

E e v := ( P £ )~ 1 Ev(P £ )~ 1 = (St)) op £ . (11) 

We assume that C s , b £ , and f £ are such that 


C e op £ = C, P b o f = b, H £ op £ = H, P g op =g, 


and that 


{f E ,v) H i/ 2 (dNns) = for every v G H^ 2 (d N tf), 


'uSo (.a N n) 


( 12 ) 

(13) 


for some coercive tensor field C G L°°(fl;Rg^ x3x3 ), b G I/ 2 (fl;R 3 ), H G L 2 (fl;R 3 ym), 
g G -ff 1 (fl;R 3 ) such that ( Eg)i$ = 0, and / G (Hgg 2 (cl^f2;R 3 ))'. From (JSJ we deduce 
that 

F £ op £ = H+ CE £ g = H+ CEg =■. F. (14) 


Remark 4.1 The required condition {Eg)n = 0 is equivalent to say that g is a 
Kirchhoff-Love displacement. Indeed, this assumption and also those on C e , B £ , H £ , 
and f £ could be relaxed. For instance, it would be enough to require that 

H £ op £ = H e , P £ g £ op £ = g £ , (15) 

for some H £ G L 2 (fl;R 3 ^), g e G -ff 1 (fl;R 3 ) which further satisfy 

g £ ^g in L 2 (fl;R 3 ), (16) 

for some g G L 2 (12; R 3 ), and 

H £ ^H, C E e g £ ^G inl 2 (Sl;R^ 3 ), (17) 

for some H,G G L 2 (fl; R 3x i). 

We note though that from G3 we have that 

E £ g £ -> C _1 G in L 2 (fi; R 3x ^), 

which, combined with CM. Korn inequality and Rellich compactness Theorem, implies 
that ( Eg)i 3 = 0, i.e., that g is a Kirchhoff-Love displacement, and that convergence 
ca actually takes place in H 1 { fl;R 3 ). 
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Remark 4.2 For the example f e considered in Remark 13.11 with f± G I/ 2 (oij t ;R 3 ) 
and // G L 2 (9jvf2|; R 3 ) which satisfy 

P e f± o F =: e/*, ° / =: /*, 

/or some G I/ 2 (w ± ;R 3 ) and fe G L 2 (8nT1 1 ', R 3 ), the rescaled surface load f defined 
by (EJ is given by 

- L /+•«*+/_/-■”*+ j w , /<•»*• 

We define the Rescaled Primal Problem as 

inf P e (v), 

where the set of rescaled admissible displacements and the rescaled energy are defined 
by 

A e := {n G R 3 ) : v = 0 on 8 d FL}, 

and 

P e {v) ■= J^CE £ vE e v - F-E e v + b-vdx-(f,v) H i, 2 i9Nny (18) 
With the assumptions m-m we have 

P{v) = eP e (v), 

where the relation between v and v is given by 00- 


Remark 4.3 In the line of Remark 1 3. SI we now make a comparison between the 
rescalings adopted for the “generalized forces” and the “standard forces”. The rescaled 
“generalized force” H contributes to the primal energy, see 00 and with the 

term 


H ■ E £ v dx, 


while the “standard forces” contribute with the terms 


(19) 


j n b-vdx+{f,v) Hl /2 {dNny 


In order to make a comparison we need to rewrite the contribution of the “standard 
forces” in a form similar to (1191) . As in Remark \3.2l given f G (Hqq 2 (8nTI-, R 3 ))' and 
b G L 2 (n e ;R 3 ), we may find H G R(div, fi) such that 


I 


b ■ v dx + (f, u) H i/ 2 (8jvfl) = I H-Evdx 


-f 


( 20 ) 


for all v £ A £ . The right hand side of © may be rewritten as 


H ■ Ev dx = / P e HP e ■ (P e )~ 1 Ev(P £ )~ 1 dx = / P e HP e ■ E e vdx, 


and the last term is exactly in the form of (1191) . Since (P e HP e )i 3 -4 0 in L 2 (f2; R 3 ^), 
while, in general , Hi 3 /0 toe deduce that the scaling of the “standard forces” is weaker 
than that applied to the “generalized forces”. 


We now change variables to the dual problem. Setting 

S :=Sop e , 

for any S £ S £ and v £ A £ we have, from m and CD. that on one hand 

[ {S - F £ ) ■ Evdx = e [ (S op* - F) ■ E e vdx 
Jne Jn 

= e [ (P £ )~ 1 (S - F)(P e y 1 ■ Evdx, 

Jo, 

— —e J div {(P e )~ 1 (S — F)(P £ )~ 1 ) ■ v dx 

+ £ ((PT\S-F)(n- 1 )n,v) H ^ dNny 

while on the other hand 


( 21 ) 


/ 


(S — F £ ) ■ Ev dx = — / div (S -F £ )-vdx + ((S - F e )n, v) ,/ 2 


Ho A (SiV^) 


/ 

JQ £ 


b e -vdx+(f,v) Hl / 2(dNQe) 


( 22 ) 


= £j^b-vdx+£(f,v) H u 2(dNeiy 

Thus from the previous two equations we find that S £ S £ if and only if 
div ((P £ )~ 1 (S — F)(P £ )~ 1 ) + 6 = 0 in L 2 (fi;R 3 ), 


((PT\S - F)(P £ )~ L )n = f in (HU a {d N (l-, R 3 ))'. 

Hence, after rescaling the admissible set S e becomes 
S e := {S £ L 2 (f2;R 3 y x i) : div ((P £ ) _1 (5 - F)(P £ )~ 1 ) + b = 0 in L 2 (fi;R 3 ) 


and ((PTHS ~ WT> = f in (^(^R 3 ))'}, 


the dual energy rewrites as 


F e *(S) := I- f C~ X S ■ Sdx, 
* Jn 


and the Rescaled Dual Problem is 


inf F £ *(S). 
ses c 


Remark 4.4 With the notation above we have 

f e *{S) = £F e *{S). 

In particular, it follows that ifT £ is the minimizer of T £ *, i.e., 

T e *(T e )= inf F £ *{S), 

SeS? 

and ifT £ is the minimizer of F e *, see ©, then 

T £ =T £ op £ . 
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Let w e := P e w s op 5 be the rescaled displacement of the solution w e of ([6]). Then 
the rescaled stress a e := a s o p e = C E e w e associated to the solution of ©• see Remark 
roi is given by 

a e =T e + CE e g = T e + C Eg. 

Remark 4.5 The rescaled dual problem coincides with the dual of the rescaled direct 
problem. 


5 Gamma-convergence of the Rescaled Dual Func¬ 
tional 

In this section, after studying the compactness of the dual problem in the weak-L 2 
topology, we identify the Gamma-limit of the sequence of dual functionals. Moreover, 
we prove the strong convergence in the L 2 topology of the minimizers. For what follows 
it is useful to notice, see (1211) and (1221) . that S G <S e if and only if 

J(S — F) ■ E e v dx = J^b ■ v dx + (/, v)^^, (23) 

for any From (1231) it easily follows that the set S e is not empty, indeed it can 

be shown that for every e > 0 there exist S s G S e such that sup e ||S s || i 2 ( f2 ) < +oo. 
This, then implies that sup e < + 00 . 

Before stating the compactness result it is convenient to set 

KLq(LI) := {v G R 1 (12;R 3 ) : (Ev)i 3 = 0, for i = 1,2,3, and » = 0on 9 dS2}, 


and 


5 := {S G L 2 {Q- R s 3 y x i) : S i3 = F i3 , for i = 1, 2, 3, and (24) 

f n ^S F) ' Ewdx = J^b ■ w dx + (f, w) H w {dNa) 

for every w G KLo(fl)}. 

Lemma 5.1 Let S e G S e be a sequence such that sup s J-£(S e ) < + 00 . Then there 
exist a subsequence, not relabeled, and an S G S such that 

S e ->• S in L 2 (n ; R 3 y x i). 

Proof. Let c > 0 be such that C ~ 1 (x)T ■ T > c|T| 2 for a.e. x G LI and for every 
symmetric matrix T. Thus 

+00 > E £ *(S £ ) > ^c||S e ||^ 2 (ri) , 

and hence sup £ ||S £ ||z. 2 (o) < +oo> which implies that there exist a subsequence, not 
relabeled, and an S G L 2 (fl; R 3 ^) 3 ) such that 

S e ->■ S in L 2 {Ll- R 3x i). 
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Let w G A'Lo(fi) and ip G ( 7 ^( 12 ; R 3 ). Define 


OX 

Va(xi,X2,X 3 ) \= W a (xi,X2,X 3 )+£ / 1 p a {xi, X2, s) ds, 

JO 

2 r 3 

{xi,X2,X 3 ) := W 3 (xi,X2,X 3 ) + £ / 1p 3 (xi,X2,s)ds. 

Jo 


V 3 


Then v e G A e , v e —» w in iL 1 (D; R 3 ) and 


P £ v £ = (P £ ) _ 1 Pv e (P £ ) _1 


Ew + 


0 0 ipi/2 \ 

0 ip 2 /2 in P 2 (D;R 3x3 ). 


\ sym ip 3 ) 

By taking S = S e and v = v e in (1231) and by passing to the limit, we find 

f (S ~ F) ■ Ew + (S - P)e 3 ■ ipdx = J b ■ w dx + (/, w) H V 2 (gNn) - 

Since w and ip are arbitrary functions, in the respective domains, we easily conclude 
that Sg5. □ 

We now identify the Gamma-limit of the dual functionals. 

Theorem 5.2 The extended functional PJ xt : P 2 (f2;R 3 ^) ->RU {+ 00 } defined by 


J "ext (5) = 


P £ *(S) ifSeSfi 
+00 otherwise, 


sequentially Gamma-converges with respect to the weak topology of L 2 ( f^R 3 ^;) to the 
functional 


Eext(S) = 


E*(S) ifSeS, 

+00 !/SgI 2 (H;R^)\ 5 , 


where 


E*{S) :=i [ C-'S-Sdx. 
2 J n 


Proof. We need to prove that: 

a) for every S G L 2 ($ 2 ; R 3 *^) and every sequence S s G L 2 (fl;R 3 X m) such that 
S e ->■ S in P 2 (D;R 3 y x i) it holds 

liminf EL t (S e ) > Pext ( 5 ); 

e 

b) for every L 2 (D; R 3 ^) there exists a sequence S £ G L 2 (f 2 ; R 3 ^) such that S e —*• S 
in L 2 (fl;R 3 ^) and 

limSUp P e £ xt (S £ ) < Pext (S). 

£ 

We start by proving a). Let S G L 2 (fi; R 3 ) and S e G Z/ 2 (f 2 ; R 3 .^) be a sequence such 
that S s ^ S in L 2 (fl;R 3 ^). We may assume that 


liminfP| xt ( 5 £ ) = limP: xt (S £ ) < + 00 . 
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Then sup e J-e X t(S e ) = sup e J>(S e ) < +oo and hence, by Lemma 15.11 it follows that 
S £ S. By a standard semicontinuity argument we have 

liminf = lim i [ C~ 1 S s ■ S s dx > i [ C ~ X S ■ Sdx = E* (S) = T ext (S). 

£ ® 2 Jn 2 Jn 

We now prove b). Let S £ L 2 (fl; R 3 ^). We may assume that T e ^t{S) < +oo. 
Thus S GS. Define Q e : S e R by 

g £ (E) := i f C _1 (E — 5) ■ (E — 5) dx. 

2 Jn 

For every e let S e be the minimizer of Q e , i.e., 


0 e (S e ) = inf £ £ (E). 


We easily check that S e is bounded in i 2 (fi; R 3 ) and hence, up to a subsequence, 

S e ^S inL 2 (fi;R 3x ^), 

for some S £ L 2 (f2;R 3 ^). Arguing as in the proof of Lemma l5.fi we also have that 
5 £ S. Taking variations we see that S e solves the following problem 


S s £ S e , 

[ C~ 1 (S e -S)-E £ dx = 0, for every E £ £ <Sg, 
Jn 


(25) 


where, 

Sq := {S £ L 2 (D;R 3 y x i) :div((P e ) -1 5(P e ) _1 ) = 0 in L 2 (D;R 3 ) 

and ((P £ ) _1 S(P £ ) _1 )n = 0 in (H^ 2 (d N V.; R 3 ))'}. 


Let 

1Z := {R £ JL(div, fl) : R ■ Ev dx = 0 for every v £ P 1 (fl;R 3 ) and v = 0 on dofl}, 

Jn 

and note that, by (1231) . 

P £ PP £ £ Sq <=> R £ H. 

Thus, let R £ 1Z and E £ = P e RP e £ So- Then E £ —¥ E in L 2 (f2;R 3 y m) for some 
E £ L 2 (fl;R 3 ^) and, arguing as in Lemma 15. II E £ <So where 

So := (5 £ L 2 (fl; R 3 y m) : Si3 = 0, for i = 1,2, 3, and 

/ S ■ Ew dx = 0 for every w £ I\Lo(fl)}. 

Jn 

Passing to the limit in (23) we then find 

j SGS, 

i J C _1 (S' — S) • Ecte = 0, for every E £ So. 

Since S, S £ S we have that S — S £ So and hence, from the above problem, we deduce 
that S = S almost everywhere in Q. As a consequence we also have that the entire 
sequence S e weakly converges to S in L 2 (fl;R 3 y ^). 
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Let B G iL(div, f2) be such that divB = b in Zr(f2; R 3 ) and Bn = —/ in 
(R^ 2 (3jvfi; R 3 ))'. Then S £ := S E - F + P £ RP e G «Sg and from (ghj we find that 

[ C ~ 1 S £ ■ (S £ - F + P £ BP £ )dx = [ C ~ 1 S-{S £ -F + P £ BP £ )dx, 

Jo. Jn 

and hence 

/ C _1 5 e • S e dx= [ C~ 1 S £ ■ (F - P e BP e ) + € _1 5 • {S e - F + P £ BP £ ) dx. 

Jo. Jn 

Denoting by B° the limit, in L 2 ($7; R 3x3 ), of P £ BP £ and passing to the limit in the 
above equation, we deduce that 

lim F £ *(S e ) = lim [ C~ 1 S £ ■ S £ dx 
e-»0 e_>0 Jn 

= [ C~ 1 S-(F-B°)+C~ 1 S-(S-F + B°)dx 

Jo. 

= [ C~ 1 S-Sdx = F*(S), (26) 

J n 

which concludes the proof. □ 


Remark 5.3 We remark that in the second part of the proof of Theorem \5.S\ we have 
indeed shown that: for every S G L 2 (fi; R 3 ^) there exists a sequence S £ G L 2 (f2; R 3 ^) 
such that S £ —¥ S in Z, 2 (f2;R 3 ) and 

limF £ ext (S £ ) = J-ext(S’). 

£ —^0 

In fact, the strong convergence of S e follows from the weak convergence and (1261) . 

Remark 5.4 In our setting, by Proposition 8.10 of f23j, sequential Gamma-convergence 
is equivalent to Gamma-convergence. 

In the next theorem we prove the strong convergence of the minimizers. 

Theorem 5.5 Let T e be the minimizer of T e * and T be the minimizer of T*. Then 

T e in L 2 (fi;R 3 y x i), 

and 

lim F £ *(T £ ) =F*(T). 

Proof. Let T £ be the minimizer of F £ *. Then by Lemma l5.ll we have that, up to a 
subsequence, T e — k T in L 2 (Q.\ R 3 ^), for some T G S. Let S G S and let S £ G S e be a 
sequence such that limsup^Q .P £ *(iS e ) < F*(S), which exists by Theorem 15.21 Since 
F £ *(T £ ) < F £ *(S £ ), by Theorem IQ we have 

F*(T) < lim inf F £ *(T £ ) < lim sup F £ * (T e ) < lim sup F £ *(S £ ) < F*(S), 

e_>0 e—f0 e->0 

which implies that T is a minimizer of T, and by taking S equal to T, that 

lim F £ *(T £ ) =T*(T). 

e —>0 
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Since J-* has a unique minimizer we have that the full sequence T s weakly converges to 
T in L 2 (f2; R 3 ^). By convexity it then follows that T e —> T in L 2 (fi; R^). Indeed, 
we have 

lim [ C _1 (T e — T)-(T e — T) dx 
e_>0 Jet 

= 2 lim (T e *(T e ) — J C _1 T e ■ T dx + E* (T)) = 0, 

from which the strong convergence follows. □ 


Remark 5.6 The rescaled stress a e = T e +C Eg associated to the solution of © , see 
Remark \4.4\ strongly converges in L 2 (fi;R 3 ^) to a ~ T + CEg. 

The next lemma, similar to a result contained in [15], allows us to characterize the 
minimizing stress tensor. 

Lemma 5.7 Let D be a bounded, open subset of R 3 with Lipschitz boundary 8D. Let 
8 dD ^ 0 be the union of a finite number of open connected sets of dD. Let 

KLo(D) := {u £ H 1 (D-, R 3 ) : (Ev)i 3 = 0, and v = 0 on 8 dD}, 


K, := {E £ L 2 (D-,M.sym) :3 z £ KLo(D) and ip £ L 2 (D;R 3 ) such that 


E = 


( Ez) a g ipp 

ip a 1p3 




and 


M = {S £ L 2 (D; R 3 y x i) : S i3 = 0, and 


/ S ■ Ezdx = 0 for every z £ I\Lo(D)}. 
J D 


Then 


K. = M . 


Proof. We first note that K, is a closed subset of L 2 [D\ R 3 yl 3 ) . Indeed, let {E-’} C 1C be 
such that E■* —¥ E in L 2 {D ; R 3 ^m). Then there exist z •' £ KLq(D) and ip-' £ L 2 (D; R 3 ) 
such that 

(■ Ez 3 ) aP (E) a ,3 ipl -¥ tpi = (E)i3, in L 2 (D), 

for some ipi £ L 2 (D). Thus to show that K, is closed it suffices to show that there 
exists a z £ I\Lo(D) such that ( E) a p = ( Ez) a g ■ But since z J £ KLo(D) we have 
that Ez 3 is a Cauchy sequence in L 2 (D; R 3 ^) and hence, from Korn’s inequality we 
deduce, in the components of D whose boundary contain part of 8dD , that z 3 —»• 3 in 
the H 1 norm, while on the other components it is z 3 minus its orthogonal projection 
on the set of infinitesimal rigid displacements which converges to some z in the H 1 
norm. Throughout D we then have ( E) a p = ( Ez) a p■ 

The proof of the lemma now follows easily. In fact, we have 1C C M ± and K. 1 - C M . 
This latter inclusion implies that C (/C x ) x . Hence 

KCM X C (£ x ) x , 

but since K. is a closed subset of L 2 (D; R 3 ^m) we have that (/C x ) x = 1C. □ 
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Theorem 5.8 The minimizer T of T* satisfies the following problem: 


T es, 

/ C _1 T ■ E dx = 0, for every E £ So, 
J n 


( 27 ) 


where 


So := {S' £ L 2 (f2;Rsy n ?) : Si 3 = 0, and 


( 28 ) 


/ 

■/n 


S • Ezdx = 0 for every z £ A'Ao(S2)}. 
Moreover, there exist a unique if £ L 2 (fl;R 3 ) and a unique u £ KLq(Q) such that 


T = ( 


( Eu) a p ipp 

Tpc, i>3 


(29) 


Proof. Problem (1271) is simply the Euler-Lagrange equation of the problem infs 6 s T* (S). 
From (EH we have that 

C _ 1 T£ (<S 0 )\ 

and hence from Lemma [5. 7 1 we deduce that there exist u £ ATLo(fi) and ip £ A 2 (12; R 3 ) 
such that 


L T = 


( Eu) a p i/j /3 
Ipc, 1p3 


Remark 5.9 The stress a = T + C Eg, limit of the stresses associated to the solutions 
of ©, see Remark roi is given by 


( (Eu + Eg) afi iff, \ 

V V’s )' 


Setting 

w := u + g £ KL g (Tl) := £ H 1 ( 12; R 3 ) : (Ev)i 3 = 0, and v = g on 9d12}, 


we may write 

a = c( {Ew)aP t' 3 1 • 

\ fpc 1p3 J 

The rescaled stresses a e = C E e w e , see Remark \4-4\ strongly converge in L 2 (12; R 3 ^) 
to a, see Remark roi thus 


( ( Ew) a p ip p \ 

V )’ 


in L 2 (f2;R 3 y x i). 


6 The bi-dimensional limit problem 

The limit problem obtained in Section[5]is defined on a three-dimensional domain. The 
aim of this Section is to show that it can be rewritten on a two-dimensional domain. 
For a given S £ S let 

n f 112 m r 1/2 

S N := / Sapdxse a ®ep, and S M := / X 3 S a p dx 3 e a ® ep. 

J - 1/2 J-1/2 
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Similarly, using the components F a p, we define F N and F M . Let 
Hq D (oj-, R 2 ) := {77 G _ff 1 (a;;R 2 ) : 77 = 0 on cfow}, 

and 

Hq ,d(w) := {v G H 2 (uo) : 77 = d n ri = 0 on 8 duj}. 

For every 2 G KLq(Q.) there exist (771,772) G H/ D {ur, R 2 ), 773 G H/ d (uj) such that 

f z a (x i,x 2 ,x 3 ) = r/ a {x i,* 2 ) - X 3 d a T) 3 (xi,X2), 

\ Z 3 (xi,X2,X 3 ) = 7 ? 3 (a:i, *2). 

A simple calculation shows that 


= {(Ep) af 3 - x 3 d a dprt 3 )e a ® e/s, 
and hence the condition, appearing in the definition of <S, see 


/ 


/ 


(S-F)-Ezdx = I b-zdx + {f,z) H ,W2 (9nQ) , 


for every 2 G KLo(ft), rewrites as 

J ( S N _ F N ) a p(E V ) a p - (S M - F M ) a pd a d pV3 dx = W n (( V i,V 2 )) + W m (7? 3 ), 


where 


and 


W ((771,772)) := 


■=\r 

Jui J -1/2 


b a dx 3 r/ a dx + (/«, ??a) H i/2 ( a N ny 


r r 1/2 

(773) := / / 63 d* 3 773 dx + (/ 3 ,773) 1/2 

Ju,J-l/2 00 


J UJ j-1/2 

ri/2 


w J-1/2 


(d N &) 


x 3 b a dx 3 d a r ) 3 dx + {f a , x 3 d a rj 3 ) H 1 / 2 (d N n)- 


Remark 6.1 If f is as in Remark \4.S[ then the work done by the loads can be written 
more explicitly, for instance 


(ftX’Va) H F 2 (g N n) ~ 
We therefore have 



r/ a dx + 



rja dx + 



dr3 r/a dx. 


5 := (S G L 2 (fi; Rsym) : 5 i3 = fi 3 , for i=l, 2 , 3 , and 


( 30 ) 



F n ) ■ Etpdx = W N (ip) for every 73 G i/o jD (w;R 2 ), 
F m ) ■ Wipdx = W m (t p) for every ip G Hq D (u}).} 


We now rewrite the functional J r * in terms of S N and S M . To do so we let 


£:={SG L 2 {Q- R^) :BA, B G L 2 (u;; R 2 *™) such that 

S a p(xi,X 2,X 3 ) = A a p(xi,X2) + l3B^(ll,*2)}- 
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Since £ is a closed subspace of Z, 2 (fl; Rg^) we have 

L 2 (n;R^) = £e£ X 

We note that E € £ x if and only if T> N = E M = Ei 3 = 0. Let II be the projection of 
Z, 2 (fl; Rs^m) to £. Then, from the relation 

[ n(S) ■ E dx = [ S ■ E dx for every E G £, 
in in 


we infer that 

n {S) a p = + i2x 3 s%,, n(5) i3 = s*. 

Hereafter we denote by 


and by 


S C := n(S') 
<S £ := n(<s) 


and 

S c 

CO 

1 

II 

and 

S c 

:=S-S C 


Lemma 6.2 With the notation just introduced we have that 

S c = £ x . 

Proof. From the definition of S c it immediately follows that S c C £ x . To prove the 
opposite inclusion first note that 

S c C 5. (31) 

Indeed, let S c G S c . Then there exists S G S such that S c = 11(5'), that is 
(, S c ) i3 = S i3 = F a , and since n{S) af> = S^+12x 3 S^ we have also that ( S C ) N = S N 
and ( S C ) AI = S M . Hence (1311) follows from the representation of S given in (1301) . 

Let E G £ x . Let S c be any element of S c . The condition E G £ x implies that 

E n = E m = S i3 = 0 and hence we have, using (1311) . that E + S x G S. Since n(E) = 0, 

n(S x ) = S c and the linearity of n, which holds because £ is a closed linear subspace, 
we have 

s = e + s c - (n(E) + n(s ,£ )) = s + s c - n(s + s £ )g5\ n(s) = s c , 

and hence £ x C S c . □ 

We may therefore write 

T*(S) = i [ C~ x S c • S c + 2C _1 S £ • S c + C~ 1 S C ■ S c dx 
Jn 

= F ± {S C ,S C ) + i [ C ~ 1 S c -S c dx, 

* Jn 

where we have set 

F ± {S c ,S c ) := [ C~ 1 S c -S c + ^C~ 1 S c -S c dx. 

JCl * 

Thus, thanks to Lemma m we have that 

inf T*(S) = inf inf F ± (S C ,S C ) + l f C -1 ^ ■ S c dx, 
se5 2 J n 
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and setting 


/ X (S £ ):= inf T ± (S C ,S c ), 
s c ec ± 

we have 

mf^(S) = g mf jC Ji ( S £ ) ) 

where we have set 

^(S' £ ):=i [ C ~ 1 S C ■ S c dx + f ± (S c ). (32) 

2 Jo 

It is possible, even for a generic elasticity tensor C, to write the function / x explicitly, 
but, as it can be seen in the next Theorem, the explicit form of / x is quite involved. 

Theorem 6.3 Let 



C aP'yS ■— ^Fafi'yS ^ 

(33) 


r 1/2 ■ - 



:= / x\Cdx 3 for i = 0,1,2, 

(34) 


J -1/2 



C := 12(C (2) - C (1) (C (0) ) _1 C (1) ), 

(35) 

and 

r := (€ (0) ) _1 + 12(€ (0) ) _1 € (1) C _1 C (1) (C (0) ) _1 , 

(36) 


C" m := -12(€ (0) )” 1 C (1) €“ 1 , 

(37) 


C m " := —C _1 C (1) (€ (0) ) _1 , C mm := C _1 . 

(38) 

For a given 

S c € S c , let A € £ x be the minimizer of F ± (S c , ■), i.e.. 



f ± (S c ) = inf F ± (S C ,S C ) = F ± {S C , A). 



s c ec ± 


Then A = C Z — S c where Z = Z + z Q e 3 , with Z := Z N + 12^3 Z M , 



( Z N :=C nn (S c ) N + C nm (S c ) M + z n , 
i Z M := C mn (S C ) N + C mrn (S C ) M + z m , 
l Zi ■■= 3 - (CZ)a), 

(39) 

and where 



z m := 

C _1 (C (1) (C (0) )“ 1 f - ff M ), z n := -(C (0) ) _1 (ff N + 12C (1) z m ), 

(40) 

with 

^a/i3 • Fi3 . 

(41) 

Moreover 

/ X (S £ ) = I [ c Z ■ Z -cr's 11 ■ S c dx + c, 

(42) 


where the constant c depends only on Fi 3 and C. 

The proof of the Theorem is given in the Appendix at the end of the paper. 
Remark 6.4 We note that: 
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1 . if C Q /3-y3 = C Q 333 = 0 , i.e., triclinic symmetry, then 


ocfi'yS 


Dt^'yS 


^ a / 333 ^ 33^5 
C3333 


2. if the material is triclinic and F a 3 = 0 then A a 3 = 0, i.e., the shear stresses are 

equal to zero. Indeed we have A a 3 = C a 3 jkZjk - 3 = 2C a 3/33-Z/33 = 2C a3 p3Zp, 

but since £p 3 = 0 it follows that zp = 3 — ( CZ ) a 3 ) = — <t^lfC a 3 ^sZ^s = 0; 

3. if C(xi, X 2 , ■) is even, for almost every ( xi,X 2 ) G oj, then C 1 - 1 ^ is null, and hence 

( Z N ■- {& 0) )-\S c ) N - 

l z M ; = t l(c( 2 ))- i (s c ) m - i(c 


4 - if Fi 3 = 0 then ff, z" and z m are null matrices; 

5. if C is independent of X 3 and Fi 3 = 0 then items @ and [^] of the present 
Remark hold and moreover 

/ x (5 ,£ ) = - [ C~ 1 S c ■ S c - C~ 1 S c ■S c dx + c. (43) 

2 in 

In fact, under these assumptions, we find C® = C, Ci 2 - ) = yjC and hence 
Z N = C“ 1 (5 £ ) jv , Z m = C from which it follows that 

Z = C-\{S c ) N + 12 x 3 {S C ) M ) = C~ 1 S" c . 

Thus from the equation of f ± given in Theorem \6.3\ it follows the representation 
of f ± given in ([43]) . The constant c, see Appendix, is equal to zero if Fi 3 = 0. 
Thus under these assumptions we have that, see (El), 

Ff(S c ):= i [ C~ 1 S c ■ S c dx. 

* Jn 

Let T c be the minimizer of Ff. i.e., 

Ff(T C )= inf Ff{S c ), 
s c es c 

and T c G S c be the minimizer of F^{T C , •), i.e., 

F ± {T C ,T C ) = inf F ± (T C ,S C ), 
s=e£ i 

then the minimizer of F* is 

T = T c + T c . 

We note that once T c is known one can determine T c directly from Theorem 16.31 
We conclude the section by noticing that the functional Ff, despite its appearance, 
it is essentially defined on ui. 


19 



7 Appendix 

This appendix is devoted to the proof of Theorem l6.3l Let S c be given and let A G 
be the minimizer of J- ± (S C , •), i.e., 

f ± (S c ) = inf F ± (S c ,S c ) = J r± (S c , A). 


Then A satisfies the following problem: 


/ 


", 1 (S C + A) • E dx = 0, for every E G C 1 


that is 


Z := C _1 (S £ + A) G C. 

Hence C Z = S c + A and since A G £ x we have that 

f (C Z) a = {S c ) i3 = F a , 

^ (C Z) N = (S C ) N , 

{ (C Z) M = ( S C ) M . 

We now show that system m delivers Z uniquely. Let 

Za •— 2^q,3, Z 3 •— Z 33: Z — Z a pCa © C/3 , 

then we havcQ 

Z = Z z (J) e 3 . 

The first equation of m rewrites as 

(CZ) i3 + (C z © es)i3 = F i3 

and by denoting, see El, 


(44) 


fcij • ^i3j3i 


it can be rewritten as 


(<z.z)i = Fi3 — (CZ) i3 . 

Since C is positive dehnite we have that c is also positive definite, and hence 

Zj = <L^(F l3 -(CZ) i3 ). (45) 

We now evaluate the in plane components of CZ. We have 

(CZ) a/3 — ^-'a/3'yS Z^S ~b ^Zaf3j3Zj — Z,yfi T Cocf3j3^ji (-^i3 (C^)^) 

— (Cck/3^5 ^Faftj3^ji ^Li3^(5) Z^fi -f Fi3- 

Setting, see El and 03, 


we have 


^Fafi'yS • — da/3'y<5 ^Fa/3j3^-ji ^-'i3 , y5 , ^a/3 ■— ^a0j3^ji F/ 3 , 


(CZ) a p — (C Z) a p + f a p. 


(46) 


4 a 0 6 := = (a ® b + b 0 a) 
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But, since Z £ C, we can write 


Z = Z N + 12 x 3 Z m , 

and, with this position, the second and third equations of flU rewrite as: 

J C (0) Z N + 12C W Z M = (S c ) N -f N , 

{ C W Z N + 12C {2) Z M = (S C ) M - ff M , ^ ’ 

where we have set, see m, 

r 1 / 2 . _ 

C 1 := / x l 3 Cdx 3 fori = 0,1,2. 

J - 1/2 

Thanks to Lemma rm below, we have 

Z N = -12(C (0) ) _1 C W Z M + (€ (0) )" 1 ((5 £ ) jv - ff^), (48) 

and 

Z M = C~ 1 ((S c ) M - ff M - € (1) (C (0) )“ 1 ((S £ ) jv - (49) 

where <C is defined by (1351) . 

Lemma 7.1 Let cc > 0 be a constant such that 

essinf <C{x)A ■ A > cc|4l| 2 , 

x£_CL 


for every symmetric matrix A G R 3x3 . 

With the notation introduced above we have 


CA ■ A = min C(A + b 0 es) • {A + b 0 es) > cc|4l | 2 

bem 3 


for every symmetric matrix A G R 2x2 . The minimum is achieved for 6 ™ m = — c - j 1 (Cj 4 )j 3 , 
and 

(C(A + 6 min © e 3 ))i 3 = 0. 

Also 

CA-A>cc\A \ 2 , 

for every symmetric matrix A G R 2x2 . 


Proof. The statements concerning C follow by an easy computation. To prove the 
statement concerning C note that 



C (B + X 3 A) ■ (B + X 3 A) dx 3 > cc 



dx3 > j^\A\ 2 , 


and since 



C (B + X 3 A) • (.B + x 3 A) dx 3 


C (0) B ■ B + 2C (1) A ■ B + C (2) A ■ A, 


we have that 

„ _ _ r 1/2 

C A ■ A = 12 min / C(B + *3^4) • (B + X 3 A) dx 3. 
SGIlsym J - 1/2 
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□ 


Thus, from (l48l) and (l49l) . and using (Ell), E3, El and El, we deduce El- 
Hence from (1391) we find Z and from (1451) we find 2 . Thus also Z is completely known 
and hence, from the relation, C Z = S c + A, also A is known in terms of S c . 

We now compute f ± (S c ). We have 

f ± (S c ) = J r± (S c , A) = E ± (S c ,CZ - S c ) 

= i [ CZ- Z-C _1 S £ -S c dx. 

2 Jn 

Let us write El as follows 

z = z min + f with zf n := -€"/(€ Z)a, f/ := cJ^F a , 


then 


CZ • Z — C(Z -f 2 0 63) • (Z -f 20 63) 

= C (Z + z min 0 e 3 ) ■ (z + z min © e 3 ) 

+ 2 C(Z + 2 mm © e 3 ) • f © e 3 + C f 0 e 3 • f 0 e 3 
= CZ-Z + € f-f, 

where to obtain the last equality we have used Lemma [7.1 1 Thus 
f ± (S c ) = i [ CZ ■ Z — C -1 ^^ • 5 ,£ + cf • fete, 

2 Jfi 

which is equivalent to El- 
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